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Convection thermo - vibrationnelle d'un systeme liquide a deux phases et stratifie 

Resume : La reponse d'un systeme liquide a deux pahses et stratifie sujet a une vibration parallele a un 
gradient de temperature impose est analysee en utilisant une methode hybride thermique Lattice- 
Boltzmann (HTLB). Les vibrations considerees correspondent aux translations sinusoidales 
d'une cavite rigide a une frequence fixe. Les couches sont couplees thermiquement et 
mecaniquement. L'interaction entre la convection thermique causee par la gravite et les 
vibrations est etudiee. La capacite des vibrations a intensifier l'ecoulement, le transfert de chaleur 
et la deformation de l'interface est etudie. Pour la gamme des parametres consideres, les resultats 
indiquent que l'effet du nombre vibratoire de Rayleigh et de la frequence vibratoire sur un 
systeme liquide a deux phases et stratifie est tres different que son effet sur un systeme liquide a 
une phase. La compraraison de ces deux systemes et discutee. 
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Abstract : The response of a two-phase stratified liquid system subject to a vibration parallel to an imposed 
temperature gradient is analyzed using a hybrid thermal lattice Boltzmann method (HTLB). The 
vibrations considered correspond to sinusoidal translations of a rigid cavity at a fixed frequency. 
The layers are thermally and mechanically coupled. Interaction between gravity-induced and 
vibration-induced thermal convection is studied. The ability of applied vibration to enhance the 
flow, heat transfer and interface distortion is investigated. For the range of conditions 
investigated, the results reveal that the effect of vibrational Rayleigh number and vibrational 
frequency on a two-phase stratified fluid system is much different than that for a single-phase 
fluid system. Comparisons of the response of a two-phase stratified fluid system with a single- 
phase fluid system are discussed. 
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1. Introduction 

Thermal convection of a single-phase fluid system in an enclosure induced simultaneously by gravity and 
vibration has received a great deal of attention in the past because of its importance in material processing or 
in heat transfer under vibrational environment [1-4] The influence of various factors, such as vibrational 
direction, vibrational Rayleigh number and frequency etc. on structure of the convective flows and the heat 
transfer in the cavity have been investigated. Particular attention has been given to mean or 
thermovibrational convection of single component fluids in plane fluid layers and cylindrical cavities, and 
includes the effects of gravity or weightlessness (See Ref.l and references therein.) 


1 



Studies of thermal convection in a two-layer fluid system with external vibration are important for 
evaluating the effect of the unavoidable vibrations of satellites and spacecraft mechanical structures on two- 
phase fluid systems. The investigation of the behavior of two-layer fluid systems subject to vibration has 
mostly dealt with isothermal systems. Kenner [5] discussed theoretically the stability of the interface 
separating two immiscible incompressible fluids of different densities and viscosities in the case of fluids 
fdling a cavity which performs horizontal harmonic oscillations. Duval et al. [6] focused on the interfacial 
dynamics of two-layer liquids under an oscillating gravitational field under isothermal conditions. They 
found that the interfacial region acts as a vortex source sheet and is susceptible to Kelvin-Helmholtz and 
Rayleigh-Taylor instabilities. The vortices produced along the interface can serve as a stirring mechanism to 
promote local mixing. Ivanova et al. [7] examined interfaces between immiscible fluids under horizontal 
vibration. 


In this paper, results of a set of simulations of vibrational convection in a non-isothermal two-layer fluid 
system are reported. The effects of external high-frequency vibration on not only the flow characteristics and 
interfacial dynamics, but also the heat transfer process are examined. 

2. Numerical Model 

The physical model is shown in Fig. 1. A rigid cavity containing two different immiscible liquids and 
subject to a periodic lateral translation, bcos(a>t)n with fixed amplitude and frequency is examined. Here b 
is the displacement amplitude, n is a unit vector parallel to the translation direction and co is the angular 
frequency of vibration. Fluid motion in the container is modeled using the two-phase thermal lattice 
Bolt z mann model (TLBM) described in Refs. 8, and 9 with a straightforward modification to account for 
motion of the cavity. Additional details concerning modeling two-phase fluids using the lattice Boltzmann 
method are found in Ref. 10. In the TLBM the fluid motion is computed using a two-particle distribution 
function approach proposed by He et al. [11, 12]. The two distribution functions satisfy the following 
equations for distribution functions £(r,e,t) and f\r,e,i) where e is a so-called microscopic velocity [8- 
10,1 1], r is the position and t is time. 

The Boltzmann equation for £(r,e,t) is [8-10] 
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The term/*) = /cpW 2 p represents the macroscopic interfacial force between the phases which is associated 
with a steep density gradient. Bulk molecular interactions are accounted for through the 
potential i// = BRT 0 p 2 % - Ap 2 , and the parameters A and B can be set so as to determine the degree of phase 
separation. For fixed values of A and B the parameter k determines the magnitude of the interfacial 
tension. The pressure and momentum are computed from distribution function %(r,e,t) as follows: 
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where the integral is over the microscopic velocity space. To calculate the density it is convenient to define 
an index function (j) (which represents, in some sense the phase fraction) such that 

4> = \f*de. (5) 

The distribution function, f* = f* ( r,e,t ) , for the index function (j) must be calculated as a function of time 
and position. This will be discussed below. The density is then calculated from the index function and is 

p(<t>) = Pi + 7 — (Ph ~ Pi)- ( 6 ) 


Here a subscript / and h denote the light and heavy fluid, respectively. Instead of setting transition values 
of the density p at T 0 = 0.883 T cr , the transition values of the index function (j) are set to be </;, = 0.0403 and 
<t> h = 0.2589 , respectively, for a given practical situation the proper density values are then determined by ( 6 ) 
after setting p h and p, to the appropriate values. The dependence of density on temperature is introduced 
through p h = p h ( T ) , and p, = p, ( 7 ) . The distribution function, f * , for the index function <f> satisfies 
[ 12 ] 
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In summary, to simulate Boussinesq-type incompressible multiphase flow, two distribution functions, c , 
are employed. The index distribution function /'’ is used to track the density distribution through 
equations (7), (5) and ( 6 ), while the auxiliary distribution function £ is used to obtain the pressure and 
momentum fields through equations (3) and (4). The lattice Bolt z mann equations are discretized and are 
coupled with a macroscopic energy equation that provides an updated temperature field which is used to 
calculate the density in the discretized body force term 6 ’(/) . 

It is convenient to describe the fluid motion in a frame of reference that moves with the translational 
displacement of the cavity. Thus, the fluid velocity, relative to this moving frame, will be zero at cavity 
boundaries. As a consequence of this change of reference frame oscillatory translations of the container 
appear as a time-dependent acceleration g(t)= g 0 -bar cos(a>t)n , where g 0 denotes the steady 
acceleration due to gravity at the earths surface. In the hybrid thermal lattice Boltzmann model (TLBM) 
g(t) appears in a body force term G(t) = p(<p)g(t) in eq. (1). 

The aspect ratio, Ar , of the cavity is set to be Ar = 2 unless otherwise stated. The upper layer fluid II 
(silicone oil, 2cSt) is assumed to be a density-linear fluid while bottom layer fluid I (water) a density- 
inversion fluid with a quadratic dependence of density on temperature. Each fluid is characterized by its 
thickness h t , thermal conductivity k t , specific heat at constant volume c vi , volumetric mass density p t , 
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kinematic viscosity v t , and thermal expansion coefficient /( . Unless otherwise mentioned, the ratio of the 
physical properties of the fluid II to fluid I have been fixed as p r =0.9, k r =0.19, c vr =0.4, v r =1.27 . 
Initially, two layers of immiscible fluids at equal heights {h x = h 2 = 0.5 H ) fill the cavity and are at rest and 
heat transfer is by conduction (a linear horizontal temperature gradient in x). The dimensionless variables 
for this problem are the Rayleigh numbers, Ra x - gA, AT 2 L? I a x v x , Ra 2 = gj3 2 AT 2 1? I a 2 v 2 (here 
AT = T h - T , where T h and T are, respectively, the hot and cold temperatures of the vertical walls), 


Pi = 


1 dp 2 


Plr 


dT 


and \ are related to the thermal expansion coefficient in the lower layer 


through = \(T -T ref ) . The Prandtl number is Pr =i /,/a 2 , and the capillary number is 
Ca = p 2 v 2 a 2 1 (7 n L , The Prandtl number is fixed at Pr = 27.6 . Another important dimensionless parameter is 
the density inversion parameter that describes the location of the density inversion temperature with respect 
to the sidewall temperatures. It is given as 
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Note that, the dimensionless temperature 6 lies in the interval 0 < 6 < 1 while -oo < I < +co . 

Two additional dimensionless parameters are necessary to characterize the vibration conditions. They are the 
vibrational Rayleigh number and dimensionless frequency and are given by Rav - ( j3 2 ba>ATL ) 2 / 2 v 2 a 2 and 
Q —cols / a 2 , respectively. The effect of translational vibration on two-phase fluid flow is examined for 
fixed Ra t = 160 , Ra 2 - 5 x 10 3 , Ca = 3.2 x 10~ 3 , I - 0.43 (thermocapillary forces are neglected). 


The local Nusselt number and the average value at vertical walls are 
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The Rayleigh numbers are kept constant at Ra t = 160 , Ra 2 - 5 x 10 3 , and Ca =3.2 x 1 0 , and / = 0.43 . 

We first calculated thermal convection in both layers due to gravity alone and then examined the effects of 
vibration for a range of vibration amplitudes and frequencies. 


3. Results and Discussion 


3.1 Convection due to gravity alone 

In the upper density-linear fluid layer, because the density gradient is negative, dpi dT < 0 , the fluid with 
lower density rises near the hot sidewall while that with higher density descends near the cold sidewall, 
which produces a clockwise convective roll (also called “regular roll”). However, in the lower layer, 
/ = 0.43 divides the fluid into two sections, the left and right sections. In the left section, the temperature is 
higher than the density inversion point and “regular roll” cell shown in Fig. 2(a) develops. In the right 
section, the temperature is lower than density inversion temperature and a second anticlockwise convective 
roll (so called “inversion roll”) cell forms. However, this flow pattern is different with that in the single- 
phase fluid with density inversion where both rolls extend from bottom to top of the liquid layer and can be 
symmetric [17]. Due to the strong viscous coupling between the top density-linear fluid layer and bottom 
density-inversion fluid layer, the convection in the top layer enhances the ’’inversion roll” cell of the right 
section and squeezes and pushes the “regular roll” cell in the left section to the left-bottom comer. A density 
inversion plane can be assumed between these two sections where / = 0.43 as shown in Fig. 2(a). The 
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temperature profiles can also reflect the effect of different vorticity convection rolls. As show in Fig. 2 (c), 
the temperature distribution in the upper layer distort seriously because a strong clockwise convection. 

3.2 Convection induced by both gravity and vibration 

Vibration in a direction parallel to the temperature gradient imposed on the two-layer fluid system produces 
the response shown in Fig. 3. The vibrational Rayleigh number and frequency for this case are Rav = 5000 
and ft = 500 , respectively. Results are presented in terms of dimensionless variables. The quantities L , 
Is / a, , AT = T h - T c are taken to scale length, time and temperature, respectively. For a high frequency 
vibration, the velocity and temperature always oscillate with a single dominant frequency as shown in Figure 
3 and the phase difference with the driving force is either close to n or n 1 2 , for velocity or temperature, 
respectively. Figure 4 displays the time histories of the average Nusselt number at the vertical cold wall. In 
general, the average Nusselt number increases from zero to a steadily oscillating periodic value. 

Figure 5 (a) - (g) shows the streamlines at different times in one oscillation period after the flow field 
reaches a steady oscillation state. The periodic flow is established quickly after about 100 periods. Within 
contrast to the case of fluid flow without vibration (see Figure 2), only one oscillating convective roll is 
generated in the entire two-layer system, as shown in Figure 5 (a) - (g). The direction and intensity of the 
convection depend on time. To further study the effect of vibration on fluid flow, it is useful to examine the 
time average values of some variables, such as velocity, temperature etc. . The mean value of the velocity is 

M=/ r J udt , (12) 


where f r is vibrational frequency f r = QHn , and the mean stream function contours are derived from the 
time averaged velocity field at steady oscillation state. 

The mean stream function at the steady oscillation state is established and shown in Fig. 5 (h). Compared 
with that of fluid flow without vibration as shown in Fig. 2 (a), the roll cell located at the left-bottom comer 
moves up to the position near the interface while another anticlockwise roll cell is generated above it. Figure 
6 (a) - (c) shows evolution of the associated temperature field. The distortion of the isotherms during the 
oscillation is small, but measurable. The mean temperature contours are plotted in Fig. 6 (d). The isotherms 
have maximum distortion near the top of the cavity and straighten out near the interface. Compared with 
Figure 2 (c), the effect of vibration is that the mean temperature field has evolved toward a conduction-like 
profile. This observation is consistent with the results in [1] and [4]. 

3.3 Effect of vibrational Rayleigh number Rav 

Figures 7 and 8 show the mean streamlines and isotherms at a fixed vibrational frequency Q = 500 for 
different vibrational Rayleigh numbers Rav. In Fig.7, two convective rolls with different senses of rotation 
near the interface and the left boundary wall grow with increased vibrational Rayleigh number Rav while the 
other two convective rolls are suppressed. The magnitude of the time-averaged velocity vector increases 
with larger values of Rav . A measure of the time-averaged velocity U is the average value of the mean 

velocity over all the nodes in the whole system 

U mg =-t J ^ + K ■ ( 13 ) 

n , 

Flere n is the total grid number, u x , u are the components of mean velocity vector li . The average velocity 
U increases with increased vibrational Rayleigh number Rav at any constant frequency Q , although the 
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rate is very small at low frequency (~!2 = 200). The increased average velocity induced by vibration 
enhances heat transfer, and results in an increased Nusselt number Nu (see Fig. 9). 

The instantaneous location of the interface also oscillates with the external vibrations, and has the 
appearance of a traveling wave. Larger values of the vibrational Rayleigh number Rav accentuate the 
distortion of the interface, as shown in Figure 10. When the vibrational Rayleigh number Rav is large 
enough, the interface can reach the top and/or bottom of cavity, and the two fluids start to mix, further 
investigation of this is beyond the scope of this paper. 

Previous studies have predicted that the thermal vibrational fluid flows of single-phase systems are very 
different from the behavior exhibited above. For a single-phase system with a vibration direction parallel to 
temperature gradient, increasing Rav will generally decrease the mean velocity and the heat transfer rate 
(Nusselt number). This is because the vibration acts as an effective gravity force that tries to reorient 
isopycnic surfaces to be perpendicular to the vibration direction [4], The mechanism for this difference is 
very interesting and is discussed below. 

For a thermal vibrational convection in an enclosure filled with single -phase fluid system under 
weightlessness ( g = 0 , Ra - 0 ), no mean motion in the fluid is expected when the vibration direction is 
perpendicular to the isotherms for high frequency vibration [4], since there is no pre-existing gravity- 
induced buoyancy flow. The situation is analogous to the unidirectional oscillatory translation of a simply 
supported rigid pendulum in a direction parallel to its length under zero gravity conditions. Flowever, in a 
two-layer fluid system at the same conditions as above, the vibration-induced phase-buoyancy force causes 
the entire system to oscillate. This is simply because the interface separating the two different density phases 
is initially flat and parallel to the vibration direction; the density gradient across the interface is 
perpendicular to the vibration direction. Thus, the interface is an isopycnic surface and will tend to rotate to 
be perpendicular to the vibration direction. This is in contrast to thermally induced density variations within 
each phase for which the local density gradient is parallel to the vibration direction. As the direction of 
translation (corresponding to the vibration) of the cavity changes from left to right, the interface will change 
shape or sloshes back and forth, sustaining an oscillatory motion of the fluid. This flow will not occur for a 
similarly oriented single-phase fluid. Thus, in the single-phase case studied by Zhao [4] with vibration 
parallel to the applied thermal gradient the mean flow approaches zero as the vibrational Rayleigh number 
Rav increases at a constant frequency Q . In contrast, for the mean flow velocity in the two-phase case 
with vibration parallel to the applied thermal gradient, the mean velocity does not approach zero but in fact 
increases with increasing values of Rav . As a result, the heat transfer rate Nu is enhanced by the vibration. 

Convection in this two-layer fluid system arises due to two sources (see Fig. 1 1); one is natural convection 
due to steady gravity oriented perpendicular to the applied temperature gradient, the second is vibration- 
induced convection which for this two-phase situation involves thermovibrational buoyancy within each 
layer and interfacial buoyancy due to the density difference between the phases. The mean flow in a two- 
layer system is enhanced by increasing the vibrational Rayleigh number Rav due to buoyant flow caused by 
the sloshing of the interface. When Rav is large enough, the mean flow induced by vibration will dominate 
the fluid flow and heat transfer in the two-layer fluid system. That is why, for the particular geometry 
studied here, a lager value of vibrational Rayleigh number intensifies the fluid flow in two-layer fluid system 
while suppresses flow in single -phase system. 

3.4 Effect of vibrational frequency D 

Figure 9 reveals the effect of vibrational frequency on the thermovibrational dynamics of a two-layer system 
on heat transfer rates. For the range of frequencies studied here 200 < Q < 1000 , and at a fixed vibrational 
Rayleigh number Rav the mean flow velocity continues to be enhanced as the frequency is increased. 
Hence, the heat transfer rate Nu is increased even though the temporal and mean distortion of the interface 
is suppressed. This is because the rapid small amplitude interfacial oscillations generate a strong mean flow. 
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The interface is a vibrational buoyancy source activated by the vibration that is roughly parallel to the 
interface. In contrast, in a single -phase fluid system, the absence of the interface allows the isotherms to 
become almost parallel to the vertical walls. In this case, the mean flow velocity and heat transfer rate are 
independent of the vibrational frequency Q [4]. 

4. Conclusions 

The effect of external vibration on the convective flow and heat transfer in a two-layer fluid system of 
immiscible liquids with density inversion was investigated using a hybrid thermal lattice Boltzmann method. 
Under the external vibration corresponding to an oscillatory uniform horizontal translation of the rigid 
cavity, oscillatory convection is established after a steady oscillation state is reached. For the range of 
conditions investigated, the results reveal that the effect of vibrational Rayleigh number Rav and vibrational 
frequency Q on a two-layer fluid system is much different than that for a single -phase fluid system when 
the vibration direction is parallel to the temperature gradient. In a two- layer system consisting of immiscible 
liquids with a density inversion, the mean flow velocities and heat transfer rate increase as Rav and Q 
increase. This is because the mean fluid motion in a two-layer fluid system subject to translational vibration 
is not zero even under weightlessness. This is in contrast to a single-phase fluid system, where the natural 
convection is suppressed and the Nusselt number decreases and tends towards unity as Rav and Q are 
increased. The interface between two-layer fluids experiences larger distortion at larger values of Rav and 
smaller values of Q . 

The effect of Rav and Q on a two-layer fluid system is much different than for a single-phase fluid system 
when the vibration direction is parallel to the temperature gradient. In a two-layer system consisting of 
immiscible liquids with a density inversion, the mean flow velocities and heat transfer rate increase as Rav 
and O increase. This is because the mean fluid motion in a two-layer fluid system subject to translational 
vibration is not zero even under weightlessness. This is in contrast to a single-phase fluid system, where the 
natural convection is suppressed and the Nusselt number decreases and tends towards unity as Rav is 
increased. The interface between two-layer fluids experiences larger distortion at larger values of Rav and 
smaller values of Q . 
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Fig. 1 Schematic of a two-layer system heated from the side walls. 
Fig.l Systeme a deux couches chauffe lateralement. 





(a) streamlines 


(b) isotherms 


Fig. 2 Buoyant thermal convection withitav = 0 . 
Fig. 2 Convection thennique pour Rav = 0 . 
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Fig. 3 Velocity ( u x ), temperature ( 6 ) and vibrational acceleration ( g v ) at node 

(25, 65) . Rav = 5 x 10 3 , Q = 500 . 

Fig. 3 Vitesse ( u x ), temperature (6 ) cl acceleration vibrationnelle ( g v ) au noeud 

(25,65). Rav = 5x 10 3 , Q = 500 . 
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Fig. 4 Time history of the average Nusselt number at the cold vertical wall. Rav = 5 x 10 3 , 

n = 500 . 

Fig. 4 Evolution du nombre de Nusselt a la paroi froide. Rav = 5 x 10 J , O = 500 . 






Fig. 5 Streamlines at different phases of one vibrational cycle after a steady oscillation is 
established. Rav = 5x10', Q = 500 . 






Fig. 5 Lignes de courant a differentes phases d’un cycle vibrationnel apres 
l’etablissement d’un regime oscillatoire. Rav = 5 x 10 3 , i2 = 500 . 



Fig. 6. Isothenns at three different phases of a cycle during a period after a steady 
oscillation is established. Rav = 5 x 10 3 , Q = 500 . (a) 0, (b) n , (c) 1 .5n , (d) time 
averaged. 

Fig. 6. Iothennes a trois phases d’un cycle durant une periode apres l’etablissement d’un 
regime oscillatoire. Rav = 5xl0 3 , 12 = 500. (a) 0, (b) n, (c) 1.5^, (d) moyenne 

temporelle. 





(a) (b) (c) (d) 

Fig. 7 Streamlines at O = 500 for different Rav . (a) 0, (b) 200, (c) 10 3 , (d) 10 4 
Fig. 7 Lignes de courant Q = 500 pour differents Rav . (a) 0, (b) 200, (c) 10 3 , (d) 10 4 
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Fig. 9 Average Nusselt number at a vertical wall Nu as a function of Rav for different 
frequencies Q . 

Fig. 9 Nombre de Nusselt moyen Nu a la paroi verticale fonction du Rav pour 
differentes frequences O . 
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Fig. 10 Dimensionless interface fluctuation 8 as a function of Rav for different 
frequencies O . 

Fig. 10 Fluctuation adimensionnelle 8 de l’interface fonction de Rav pour differentes 
frequences Q . 




(a) thermal vibrational (b) steady gravity (c) vibration induced 

convection induced convection convection 


Fig. 1 1 Decomposition of the thennal vibrational convection mechanism in the presence 
of a deformable interface between two immiscible liquids. 

Fig. 1 1 Decomposition du mecanisme thermo-vibrationnel en presence d’interface 
deformable entre deux liquides immiscibles. 






